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Abstract: The surface wind speed (SWS) is affected by both large-scale circulation and land use and cover 
change (LUCC). In China, most studies have considered the effect of large-scale circulation rather than 
LUCC on SWS. In this study, we evaluated the effects of LUCC on the SWS decrease during 1979-2015 
over China using the observation minus reanalysis (OMR) method. There were two key findings: (1) 
Observed wind speed declined significantly at a rate of 0.0112 m/(s-a), whereas ERA-Interim, which can 
only capture the inter-annual variation of observed data, indicated a gentle downward trend. The effects 
of LUCC on SWS were distinct and caused a decrease of 0.0124 m/(s-a) in SWS; (2) Due to variations in 
the characteristics of land use types across different regions, the influence of LUCC on SWS also varied. 
The observed wind speed showed a rapid decline over cultivated land in Northwest China, as well as a 
decrease in China's northeastern and eastern plain regions due to the urbanization. However, in the 
Tibetan Plateau, the impact of LUCC on wind speed was only slight and can thus be ignored. 


Keywords: surface wind speed (SWS); land use and cover change (LUCC); observation minus reanalysis (OMR); 
normalized difference vegetation index (NDVI); China 


Citation: LI Yupeng, CHEN Yaning, LI Zhi. 2019. Effects of land use and cover change on surface wind speed in China. 
Journal of Arid Land, 11(3): 345-356. https://doi.org/10.1007/s40333-019-0095-5 


1 Introduction 


Surface wind speed (SWS) partially governs the transfer of water, energy and momentum 
between the land surface and the lower atmosphere (Azorin-Molina et al., 2014). Understanding 
the evolution of long-term SWS is of great significance for evaluating coastal erosion, surface 
energy balance, hydrological cycle, pollutant dispersion and wind energy potential (Viles and 
Goudie, 2003; Pryor et al., 2006; Roderick et al., 2007; Pineda-Martinez et al., 2011). Recently, to 
explain the observed decrease in pan-evaporation in Australia, Roderick et al. (2007) proposed the 
concept of "wind stilling". A meta-analysis of 148 regions was conducted by McVicar et al. 
(2012), which indicated that "wind stilling" widespread in the tropics and mid-latitudes of both 
hemispheres. Specifically, in mid-latitudes, the SWS has declined by between 0.004 and 0.017 
m/(s-a) over the past 30-50 a (Roderick et al., 2007). A large number of regional studies have also 
found that the prevalence of wind decline, such as the contiguous United States (1973-2000) 
(Pryor et al., 2009), the Czech Republic (1961-2005) (Brázdil et al., 2009), Switzerland (1981- 
1998) (Weber and Furger, 2001), Portugal and Spain (1961-2011) (Azorin-Molina et al., 2014), 
Germany (2007-2010) (Kaiser-Weiss et al., 2015), Australia (1975-2006; 1975-2004; 1975- 
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2006) (McVicar et al., 2008; Roderick et al., 2007; Troccoli et al., 2012), the west coast of 
Canada (begin in the late 1940s or the 1950s and run through to the early to mid-1990s) (Tuller, 
2004), South Korea (mid-1950s—2003) (Kim and Paik, 2015), and Turkey (1975-2006) 
(Dadaser-Celik and Cengiz, 2014). 

In China, changes in SWS have also attracted attention in recent decades (Xu et al. 2006; Jiang 
et al., 2009; Fu et al., 2010; Guo et al., 2011; Yang et al., 2012; Lin et al., 2013; You et al., 2014; 
Wu et al., 2016; Zha et al., 2016a; Zha et al., 2016b; Wu et al., 2017). For example, Xu et al. 
(2006) reported that the annual mean wind speed over China decreased steadily by 28% and the 
prevalence of windy days (daily mean wind speed >5 m/s) decreased by 58% between 1969 and 
2000. Jiang et al. (2009) indicated that the regions with declining trends match those with 
relatively strong observed winds, while the regions without significant declining trends match 
those with light observed winds. Lin et al. (2013) further proposed that SWS at high altitudes is 
more sensitive to climate change. 

The reason for the decrease in SWS is unclear, mainly because the changes are induced by a 
combination of anthropogenic activities and climate shifts. Anthropogenic activities primarily 
lead to land use and cover change (LUCC), while changes in the natural climate are associated 
with alterations in large-scale circulation. Thus, studies of the long-term changes and causes of 
SWS need to take into account the effects of large-scale circulations and LUCC activities as well. 

Some studies indicate that recent declines in SWS are due to changes in large-scale circulation 
patterns and that shifts in atmospheric circulation can explain 10%—50% of the decline in wind 
speeds (Vautard et al., 2010). For example, Yu et al. (2015) discovered that changes in SWS in the 
United States are closely related to the North Atlantic Oscillation in summer and the El 
Nifio/Southern Oscillation. Dadaser-Celik and Cengiz (2014) suggested that large-scale 
circulation and temperature fluctuations are the main causes of wind speed changes in Turkey. 
Across China, the decline in SWS has been attributed to a north-south warming gradient in 
winter, and to sunlight dimming caused by air pollution lingering over central areas in summer 
(Xu et al., 2006; Guo et al., 2011). The Intergovernmental Panel on Climate Change has reported 
that the poleward movement of cyclone activity is the main reason for the decline in mid-latitude 
wind speed as well as the rise in high-latitude wind speed. 

At the same time, a significant decline in SWS has been reported in the tropics and mid-latitudes 
of both hemispheres, whereas several coastal studies have reported an increasing SWS (McVicar et 
al., 2012). These reports are in agreement with observed increases in oceanic SWS measured by 
both remote sensing and in situ systems (McVicar et al., 2012; Ma et al., 2016). Therefore, some 
investigators suggested that the SWS slowdown could be caused by increased land surface 
roughness induced by LUCC (Wu et al., 2016, 2017; Zha et al., 2016a, b; Li et al., 2017). 
Mesoscale model simulations do indicate that an increase in surface roughness could explain 25%- 
60% of the observed wind stilling (Vautard et al., 2010). LUCC, including urbanization, forestation 
and decreases in pasture land area, are probable causes of the increase in surface roughness (Wever, 
2012). 

In China, under the influence of LUCC, the mean wind energy has weakened by 3.84 W/m per 
decade during 1960-1999 (Li et al., 2008) and the SWS has decreased by 0.12 m/s per decade 
during 1980-2011 (Zha et al., 2016b). Furthermore, Zha et al. (2016b) suggested that a decrease 
of 0.1 m/s in SWS could be induced by a 10% rise in the urbanization rate. By applying the 
frictional wind model, which can be used to estimate the relative effects of the pressure-gradient 
force and drag on SWS changes, Wu et al. (2016) also arrived at the same conclusion. 

LUCC is the most important and direct factor affecting SWS (Wu et al., 2016; Zha et al., 
2016a). LUCC occurs due to human impacts such as urbanization, farmland reclamation and 
vegetation destruction under the background of global climate change. When dealing with the 
impact of LUCC on SWS, much of the previous research has focused on urban land, while other 
land use types have received considerably less attention. This is unfortunate, because factors such 
as afforestation, grassland degradation, land degradation and irrigation can also change the 
surface heat flux, and thus indirectly impact on atmospheric circulation (Wever, 2012). Although 
the effects of other land use types are difficult to separate from the increased levels of greenhouse 
gases, ignoring their presence will not further our understanding of how LUCC impacts on SWS. 
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The two most important anthropogenic activities that impact on SWS changes are the increase in 
greenhouse gases and changes in land use (e.g., urbanization and agriculture), although the impacts 
of these activities are difficult to separate. In previous studies, the influence of LUCC on SWS has 
mainly been focused on urbanization, with comparisons made with rural stations. "Urban" and 
"rural" areas are generally defined according to population density (Xu et al., 2006) or light data 
(Li et al., 2017). However, this approach has a number of problems: (1) Most stations are located 
near cities, whereas only a few stations are located in mountain regions; (2) The division according 
to population has a certain degree of subjectivity and does not take into account the level of 
economic development in different regions; and (3) Although there is some objectivity in the 
relationship of lighting data with population data, the two methods are static and do not enable a 
dynamic analysis. With the development of economy, some rural stations will also be affected by 
urbanization. 

Kalnay and Cai (2003) showed that only relying on the data comparisons between the urban 
and rural stations will underestimate the role of LUCC. In response to this dilemma, they 
proposed a way to estimate the influence of land use changes over climate variables by comparing 
observational and reanalysis data using the observation minus reanalysis (OMR) method. The 
theoretical basis of the method is that the reanalysis dataset contains elements pertaining to the 
influence of large-scale climate changes caused by greenhouse gases, and while the reanalysis 
dataset is less sensitive to regional surface processes associated with different land use types, the 
surface observation data does contain all of the radiative forcing. Therefore, the difference 
between surface observations and reanalysis data can be partly attributed to surface factors such 
as urbanization and agriculture. In addition, the OMR method allows for SWS trend induced by 
natural climate variability to be cancelled (i.e., temporary changes in circulation) because they are 
present in both the reanalysis dataset and observations. The OMR approach is a useful technique 
for studying the impact of anthropogenic LUCC on climate and has been used in numerous 
applications (Wu et al., 2016, 2017; Zha et al., 2016a, b). 

In China, previous studies of the impact of land use on wind speed have mainly dealt with the 
densely populated eastern region. It is unclear whether SWS changes in more sparsely populated 
regions are also heavily influenced by urbanization. The objective of this study is to assess the 
sensitivity of surface wind speed to the influence of LUCC in China using the difference between 
the observation and reanalysis data. 


2 Dataset and methods 


2.1 Study area 


The Chinese mainland was selected as the study area (18°—53°N, 73°-135°E). The landscapes of 
the area vary significantly across its vast width. In the east, along the shores of the Yellow Sea 
and the East China Sea, there are extensive and densely populated alluvial plains, while on the 
edges of the Inner Mongolian Plateau in the north, broad grasslands predominate. Southern China 
is dominated by hills and low mountain ranges, while the central-eastern region hosts the deltas of 
China's two major rivers, the Yellow River and the Yangtze River. China is the world leader in 
wind power generation, with the largest installed capacity of any nation and continued rapid 
growth in new wind facilities. With a large land mass and long coastline, China has relatively 
abundant wind resources. According to estimates made by the China Meteorological 
Administration, based on the observations with the relatively low height at 10 m above ground, 
the theoretically exploitable wind resource represents a potential power generation capacity of 
4350 GW and a technically exploitable wind resource of 297 GW (Xia and Song, 2009). 


2.2 Dataset 


2.2.1 Reanalysis data 


Daily wind data at 10 m above ground was obtained from the latest European Centre for 
Medium-Range Weather Forecasts interim reanalysis (ERA-Interim) dataset, covering the period 
1979-2015. This dataset (Dee et al., 2011) has a resolution of 0.125°x0.125° and is used to 
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represent the long-term effects of climate change on SWS. The data assimilation system used to 
produce ERA-Interim, which is superior in quality to ERA-40, is based on the 2006 release of the 
integrated forecast system, IFS-Cy31r2 (Berrisford et al., 2011). The land surface observations 
assimilated in ERA-Interim do not include surface wind observations over land (Dee et al., 2011). 
Therefore, while ERA-Interim reanalysis dataset do include large-scale natural climate change 
signals, they do not indicate the impact of local LUCC on SWS. 

The supplementary reanalysis data used in this research were daily mean 10-m wind speed 
derived from the NCEP (the United States National Centers for Environmental Prediction)/NCAR 
(the National Center for Atmospheric Research) Reanalysis dataset (NNR), which we obtained 
from _https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.derived.surface.html. The 
NNR dataset were derived using a gaussian grid at 2.5° resolution and had the same temporal 
coverage as the observations for the period 1979-2015. Furthermore, the NNR does not 
assimilate surface observations directly, but instead estimates SWS from rawinsonde- and 
satellite-derived data. This means that the NNR estimates are also insensitive to local LUCC. 
2.2.2 Ground observation data 
The daily wind speed data from 1979 to 2015 for 602 observation stations were obtained from the 
China Meteorological Administration and compared with SWS from the ERA-Interim dataset. 
The stations used in this research were selected according to the following criteria: (1) the data 
must cover the period from 1979 to 2015; (2) any missing data must not exceed 1% of the entire 
study period; and (3) the correlation coefficient for the relationship between station data and 
ERA-Interim dataset SWS must pass a significance test and have a confidence level of 0.01. The 
linear interpolation is used to remove the stations with missing values. We applied a linear 
interpolation for stations with missing values. According to these criteria, 526 of the 602 stations 
were selected, the distribution of which is shown in Figure 1. As can be seen, the stations were 
mainly located in the eastern flat areas of China, with correlation coefficients for the relationship 
between station and ERA-interim dataset higher than 0.6. The number of stations in the complex 
terrain of western China was significantly less and the correlation coefficient was lower, 
especially in the Qaidam and Junggar basins. For all selected stations, the average value of the 
correlation coefficient was 0.44, but the correlation was lower in western China, mainly due to the 
low data density and proximity of mountains. 
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Fig. 1 Correlation of surface wind speed (SWS) between ERA-Interim dataset and surface observations. Data 
were not available for the Taiwan region. 


2.2.3 Normalized difference vegetation index (NDVI) 

The NDVI can be used as a proxy for LUCC and can also quantitatively establish the relationship 
with wind speed. The latest version of the NDVI3g.v1 product (Tucker et al., 2005), covering 
July 1981 to December 2015, was derived from data gathered by the national oceanic and 
atmospheric administration advanced very high resolution radiometer (NOAA AVHRR). NDVI 
monthly data at 8-km spatial resolution were generated from already processed 15-day NDVI 
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values. The processed data were corrected for calibration, view geometry, volcanic aerosols and 
other effects not related to vegetation change, using the maximum value composites method. A 
calibration was also applied to the original data to minimize the effects of sensor degradation 
based on invariant desert targets (Tucker et al., 2005). This dataset is considered to be the most 
accurate long-term AVHRR data record, which makes it appropriate for long-term studies of land 
surface trends in vegetation and seasonality, and for coupling climate variability and vegetation. 
The available NDVI data were averaged over nine cells surrounding each station to calculate the 
linear trends from the annually averaged NDVI. 


2.2.4 Land use data 


The land use dataset, on a scale of 1:100,000 and covering the period from the late 1980s to 2015, 
was downloaded from the Resource and Environment Data Cloud Platform (China). It was 
generated using the human-computer interactive interpretation method, which uses information 
regarding remotely sensed land-use changes to interpret Landsat Thematic Mapper (TM) digital 
images. In addition, the dataset was also confirmed by field investigation records and photographs 
(Liu et al., 2003, 2014), which contained six classes and 25 sub-classes of land use, including 
cultivated land, woodland, grassland, and urban land. The accuracy of the six classes of land use 
and the overall accuracy of the 25 sub-classes was above 94.3% and 91.2%, respectively (Liu et 
al., 2003, 2014). 


2.3 Methodology 


The ERA-Interim dataset were interpolated to the 526 selected meteorological stations using the 
bilinear interpolation method, which is suitable for converting one grid forecast field to another or 
to a discrete observation field. A non-parametric Mann-Kendall test was also used to calculate 
trends in NDVI, ERA-Interim and observation wind speeds. A trend was considered significant 
when a confidence level of 95% was achieved. At the same time, daily data was integrated into 
seasonal and annual data and the season was divided as: spring, March to May; summer, June to 
August; autumn, September to November; and winter, December to February of the next year. 

To test the influence of different land types on SWS, we first calculated the area of different 
land-use types in the two periods (1980s and 2010s) within a 10 km range of the station, and then 
calculated the rate of change of the land-use type. Then, we compared the mean value of the SWS 
trend between areas that experienced increases and stations that experienced decreases for given 
land-use types, using a paired-samples t-test. Due to the small overall proportion of water area 
and unused land, and given that the stations were rarely located on these two land-use types, we 
mainly considered the influence of changes on near-SWS in cultivated land, grassland, water 
bodies and urban land. 

A land-use dataset during 1980s—2015 was used to calculate the change in the area of different 
land-use types and the relative changes in the range of -10%—10% around a set value of 0. As can 
be seen from Figure 2, the amount of cultivated land in northern China increased, especially in 
the northern slopes of the Tianshan Mountains and the Northeast Plain, while the amount of 
cultivated land in southern China decreased (Fig. 2a). The amount of forest area was reduced in 
the key forested areas, such as northeastern China, because of logging and the expansion of 
farming areas. In the Loess Plateau and hilly areas in the south, forest land increased due to the 
implementation of Grain for Green Project in 2000 to encourage the return of farmland to forests 
(Fig. 2b). Throughout the entire country, there was a general reduction in grasslands, especially at 
the edge of the Tarim Basin, in the Northeast Plain, and on the southeast coast (Fig. 2c). 

In recent decades, the speed of urbanization accelerated and the urban land area also increased, 
especially in the heavily economically developed Beijing-Tianjin, Yangtze River Delta, and Pearl 
River Delta regions (Fig. 2d). In the northwest region, there was a decrease in grassland and an 
increase in cultivated land, whereas in the northeast, the trend was toward increased urbanization 
of land and a reduction in forestland and grassland areas. In the eastern plains, there was a rapid 
expansion in urban land, but there was little change in land-use in the Tibet Plateau. 

Based on the above characteristics, we divided the research region into four sub-regions 
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(Northwest China (I), Northeast China (II), Tibetan Plateau (III) and Eastern Plain (IV)) and 
investigated the influence of different land-use types on the changes in SWS in the different 
regions. We compared the mean trend of SWS in areas with increases and decreases in the same 
land-use type. This standard did not apply to the urban land in China because the urban land area 
has increased across most of the country. Previous studies have classified cities as either "large" 
or "small" according to population size. However, classification according to population has a 
certain subjectivity and does not take into account the level of economic development in different 
regions. Moreover, this approach is static and does not enable a dynamic analysis. In the present 
study, we only considered changes in the surrounding environment of each chosen station. We 
compared the different rates of urbanization only if the stations displayed a significant difference 
in SWS over time. This enabled a comparison of the trend in mean SWS between urban land 
areas that had expanded by greater than 10% and urban land areas that had expanded by less than 
10%. 
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Fig. 2 Distribution of national land-use changes from the late 1980s to 2015 in the four geographical units. I, 
Northwest; II, Northeast; III, Tibetan Plateau; IV, Eastern Plain. Data were not available for the Taiwan region. 
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3 Results 


3.1 Spatial-temporal variations of SWS in the observation and ERA-Interim dataset 


Figure 3 compares 10-m monthly SWS changes from observations and ERA-Interim dataset. The 
ERA-Interim dataset very well captured the seasonal and inter-annual fluctuations of SWS. The 
correlation coefficient for the relationship of the ERA-Interim dataset with observations was 0.64, 
giving a significant at the 0.99 level. The rate of decrease in SWS at the stations reached 0.0112 
m/(s-a) between 1979 and 2015. The wind speed across all seasons was significantly decreased, 
especially in spring, when the decrease in SWS was 0.0147 m/(s-a) (Table 1). However, the 10-m 
SWS of ERA-Interim dataset showed only a slight decrease of 0.0004 m/(s-a) (NNR: —0.0069 
m/(s-a)), which was non-significant at the 0.01 level. 

We also compared the mean SWS for different periods. We found that the difference between 
the ground observation data and ERA-Interim dataset increased from 0.50 m/s in the 1980s to 
0.79 m/s in the 2010s. This difference could be attributed to urbanization and other surface 
changes, because ERA-Interim dataset only reflects the influence of large-scale circulation 
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patterns induced by greenhouse gases on SWS, and is insensitive to the local and regional surface 
processes associated with different land-use types. This change corresponded to China's reform 
and opening-up period, when economic development and urbanization accelerated and the human 
impact on the land surface grew rapidly. 
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Table 1 Linear trend of surface wind speed (SWS) in China during 1979-2015 
Linear trend of SWS (m/(s-a)) 
SWS 
Annual Spring Summer Autumn Winter 
Observation —0.0112™ —0.0147™ —0.0096™ —0.0095"" —0.0096" 
ERA-Interim —0.0004 0.0200 ~0.0038" —0.0013 0.0025" 


Note: *, 0.05 level of significance; **, 0.01 level of significance. 


Figure 4 shows the spatial distribution of SWS values and linear trends from the observations 
and ERA-Interim dataset. Several conclusions could be drawn from the station observation data. 
Higher SWS values occurred in the north of the Yellow River, where the SWS was higher than 2.5 
m/s, while the SWS was generally lower than 2.5 m/s in the south of the Yellow River, especially 
in the Sichuan Basin (Fig. 4a). This finding was also reported by Jiang et al. (2009) and Zha et al. 
(2016b). The coastal SWS was greater than the inland SWS, which may be attributed to stations 
near the coast having less resistance and the general increase in wind speeds over the sea in recent 
years (McVicar et al., 2012). Due to the surface resistance being relatively insignificant, the SWS 
at high-altitude stations was larger than at low-altitude stations. Although the distribution of the 
ERA-Interim dataset wind speeds were similar to the observed maximum and minimum regions 
(Fig. 4b), the ERA-Interim dataset wind speeds were distinctly higher than the observations for 
almost all of China. In addition, the observed SWS showed a noticeable decline in most parts of 
the country. The rate of decrease in SWS reached 0.0127 m/(s-a) in 402 of the 526 stations that 
showed a wind speed decline (Fig. 4c). Moreover, there were significant decreases at 75% of the 
stations, and the regions with declining trends matched the areas with relatively strong observed 
winds. However, the 10-m wind speed of ERA-Interim dataset decreased by only 0.0003 m/(s-a), 
with the wind speed at 262 stations showing a rising trend (Fig. 4d). 

From the above analysis, it was apparent that large-scale circulation could influence the 
near-SWS but was not the main reason for the decreases in SWS because ERA-Interim dataset 
can only capture the seasonal and inter-annual fluctuations of observations and showed indistinct 
decreases. Given that SWS is affected by large-scale circulation and LUCC, we further 
investigated the influence of land cover change on SWS. 
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Fig. 4 Spatial distribution of SWS according to observations (a) and ERA-Interim dataset averaged from 1979 
to 2015 (b), and their linear trends (c, d). Data were not available for the Taiwan region. 


3.2 Influence of LUCC on SWS 


The OMR method can effectively separate the impact of LUCC on SWS from other influencing 
factors. As shown in Figure 5a, the majority of stations displayed a downward trend, especially in 
the Northeast Plain and coastal regions. This result was also confirmed by the NNR dataset. 
Strangely, the declining trend reversed to an increasing trend after the start of the 21* century 
(Fig. 5b). This phenomenon was also reported by Xu et al. (2006) and Zha et al. (2016b). Because 
changes in SWS are related to LUCC, we also calculated the NDVI value for the area surrounding 
the stations. The results showed that vegetation status had significantly improved in recent 
decades, with 85% of stations showing a rising trend in NDVI values that was closely linked to 
SWS (Fig. 6a). The relationship had a correlation coefficient of —0.743, with a significant at 0.99 
level (Fig. 6b). 

The results, as shown in Table 2, indicated significant differences in the mean SWS trend 
across all land use types except for urban land in the northwest region. The SWS declined much 
more significantly in areas where the area of cultivated land increased, largely due to the decrease 
in woodland and grassland (0. 0113 m/(s-a)), rather than in locations where the area of cultivated 
land decreased (0.0071 m/(s-a)). We therefore inferred that the increase in the area of cultivated 
land is an important reason for the decrease in SWS in Northwest China. 

There was a significant difference in the mean SWS trend over urban land in the northeast and 
eastern plain regions, where the urban land area increased by 23% and 51%, respectively. The 
SWS decline in the northeast region was greater than that in the eastern plain region, which was 
consistent with the spatial distribution of SWS changes. However, this might also be related to the 
scope of the eastern plain region selected for the study, with many upward wind speed stations 
located in the Sichuan Basin. 

In the Tibetan Plateau, there was almost no change in the areas of cultivated land, woodland, 
and grassland and there was also no significant difference in the mean SWS trend across these 
land use types. Although there was a slight increase in the urban land area, it had no real influence 
on SWS. Therefore, the effects of LUCC on SWS can be ignored in the Tibetan Plateau. 
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Fig. 5 Spatial distribution of the linear trend in SWSD (difference of SWS between the ERA-Interim dataset 
and observations (a)) and temporal variability of SWSD (b). Data were not available for the Taiwan region. 


60°E 80°E 100°E 120°E 140°E 2.8 7 SWS r 0.44 
0 500 1000km (b) —— NDVI 
40°N! 2.6 0.42 
= 
2.44 +0.40 O 
30°N} Z 
weie! 
Legend Y a TS 
Soni Laimer TEAN 2.24 0.38 
20°N NDVI trend a $ 
+ -0.023~-0.010 Porth oor, a NY 
+ -0.010-0.000 + 0. a pa; ec 
«_0,000-0.010 :gote-g07% 9 210) cee eee 0.36 


Approved No. GS(2019)2461 1975 1980 1985 1990 1995 2000 2005 2010 2015 


Fig.6 Annual trends (per decade) of average NDVI values over a 24 kmx24 km area surrounding the stations as 
obtained from the GIMMS AVHRR product (a) and temporal dynamics of annual average SWS and NDVI (b). 
Data were not available for the Taiwan region. 


Table 2 Linear trends and areal changes for different land use types for the period 1979-2015 


p Linear trend (—) Linear trend (+) 1980s 2015 Variation 
Region Land type 
(m/(s-a)) (x10* km?) (%) 

Cultivated land —0.0071* —0.0113* 6.90 9.19 33.20 

Northwest Woodland —0.0065* —0.0150* 4.13 4.08 -1.17 

China Grassland —0.0125" —0.0002* 44.62 42.58 —4.57 

Urban land —0.0101 —0.0080 0.44 0.83 89.98 

Cultivated land —0.0218 -0.0213 36.59 40.34 10.24 

Northeast Woodland —0.0222 —0.0196 52.90 51.08 -3.44 

China Grassland -0.0219 -0.0203 57.26 54.58 —4.68 

Urban —0.0196* —0.0249* 3.21 3.96 23.41 

Cultivated land —0.0087 —0.0047 128.76 124.05 -3.66 

Eastern Woodland —0.0086 —0.0076 139.00 140.27 0.19 

Plain Grassland —0.0078 —0.0102 51.58 50.47 -2.15 

Urban —0.0064* —0.0112* 10.63 16.04 50.83 

Cultivated land —0.0161 —0.0104 1.84 1.86 1.24 

Tibetan Woodland —0.0143 —0.0238 27.42 27.36 -0.23 

Plateau Grassland -0.0153 —0.0144 150.72 150.81 0.06 

Urban -0.0151 NA 0.11 0.23 105.08 


Note: —, stations with the average areal decreases; +, stations with the average areal increases; *, 0.05 level of significance; NA, not 


available. 
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4 Discussion 


Urbanization, as the most typical human reconstruction of the land surface, has received much 
attention in relation to its influence on SWS change. However, ignoring other land use types will 
underestimate the effects of LUCC on SWS. Based on the characteristics of LUCC in China, we 
divided the country into four sub-regions and considered each in relation to the impact of LUCC 
on SWS. 

The increase in cultivated land was the main cause of the decrease in SWS in Northwest China. 
During 1979-2015, the area of cultivated land in the region increased by 2.29104 km? (Liu et al., 
2014). Liu et al. (2014) reported the characteristics of cultivated land in China as "the total basic 
balance, with southern areas decreasing and northern areas increasing, and the gravity of new 
arable land gradually moving from northeast to northwest". Han et al. (2016) investigated the 
relationship between the change in SWS and area of cultivated land within a 4-km radius of 135 
meteorological stations. In Northwest China, agricultural development could explain 37.8%- 
54.2% of the observed wind stilling in the growing season, and 42.4%-—49.0% in the non-growing 
season. Han et al. (2016) summarized how agricultural development can affect SWS. They 
reported that enhanced irrigation and fertilization are conducive to the growth of vegetation, thus, 
increasing the local roughness length. With an increase in crop cultivation, the area of shelterbelt 
was also increased to improve the farmland microclimate and ensure a high and stable crop yield. 
For example, in Xinjiang, the proportional area of farmland shelterbelts increased from 2.2% in 
1977 to 5.0% in 1989 and then to 9.5% in 2008 (Zheng et al., 2013). The development of irrigated 
areas in large oasis regions may be associated with local circulations that lead to an "oasis effect", 
especially in the growing season, which can reduce the near-SWS (Ozdogan et al., 2004). 

Urbanization, which can explain the decline in SWS in northeastern and eastern China, has 
attracted much attention in relation to the influence of wind speed. For example, the SWS at 
urban stations has been shown to be lower than that at suburban stations in Beijing (Zhen et al., 
2011). Using the weather research and forecasting (WRF) model, Zhang et al. (2010) simulated 
how urbanization can induce SWS decreases of up to 50% in the Yangtze River Delta. Similarly, 
Zha et al. (2016b) showed that China's urbanization rate, which is increasing by 9.36% every 10 
years, can lead to reductions in wind speed of 0.11 m/s. 

Urbanization has many influences on SWS. Urbanization forms an urban canopy structure, 
which changes the local land cover and underlying surface radiation, thermal and dynamic 
characteristics, thus, making urban land cover very different from the surrounding areas. This 
difference affects the characteristics of the various boundary layers as well as the energy 
exchange between the land surface and the atmosphere (Zhang et al., 2010). Urbanization also 
increases the roughness of the ground, which reduces the wind speed to a greater extent than the 
surrounding suburbs. 

Over the period 1979-2015, LUCC can essentially be ignored in the Tibetan Plateau region, 
but the wind speed still showed a downward trend, which may be related to large-scale 
atmospheric changes. You et al. (2014) found that SWS in the Tibetan Plateau decreased at a rate 
of 0.024 m/(s-a) from 1980 to 2005, with the pressure gradient between high- and low- latitude 
bands changing the regional atmospheric circulation and then influencing the variations in wind 
speed. 

Additionally, reductions in snow cover and afforestation also exert an influence on the 
near-SWS, because reduced snow cover increases the surface roughness. Yet another factor 
contributing to the decrease in wind speed is the post-2000 implementation of Grain for Green 
project that encouraged the return of farmland to forest, which has resulted in an increase in the 
forest area of 2.37x10° hm? (Liu et al., 2003). 


5 Conclusions 


In this study, the influence of LUCC on SWS in China was analyzed using the OMR method for 
the period of 1979-2015. Observed SWS decreased significantly, with the largest downward trend 
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(0.0147 m/(s-a)) occurring in spring. However, the changes in wind speed in ERA-Interim dataset 
showed a relatively gentle and insignificant decline of only 0.0004 m/(s-a). The ERA-Interim 
dataset captured the seasonal and inter-annual fluctuations of the observed time series, but not the 
main reason for the decline in SWS. LUCC caused near-SWS decreases at a rate of 0.0124 
m/(s-a). The correlation coefficient for the relationship between NDVI and near-SWS was —0.733, 
Le., a significant negative correlation, which further proved that LUCC is the main reason for the 
SWS decline in China. 

In relation to the characteristics of LUCC in China, we considered the effects of LUCC on 
SWS in different sub-regions. The main reason for the decline in SWS in the northwest was the 
expansion of cultivated land, while SWS decreases in Northeast China and Eastern Plain were the 
result of urbanization. In the Tibetan Plateau, the effects of LUCC on near-SWS were so minor 
that they could be neglected. 
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